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The nucleation and growth kinetics of highly ordered nanocrystal superlattices during the evap- 
oration of nanocrystal colloidal droplets was elucidated by in situ time resolved small-angle x-ray 
scattering. We demonstrated for the first time that evaporation kinetics can affect the dimensional- 
ity of the superlattices. The formation of two-dimensional nanocrystal superlattices at the liquid-air 
interface of the droplet has an exponential growth kinetics that originates from interface "crushing". 

PACS numbers: 78.67.Bf, 61.46.+W, 81.16.Dn 



The microscopic mechanism for self-organization of 
ligand-stablized nanocrystals has been well established 
[lj, |2| • The interparticle van der Waals interaction pro- 
vides the attractive force to induce the self-assembly, 
whereas the steric interaction due to the ligand inter- 
digitation provides the balancing force to create sta- 
ble structures. However, the macroscopic patterns cre- 
ated by such interactions can vary dramatically un- 
der different experimental conditions. Even for a sim- 
ple case of evaporating a colloidal droplet of nanocrys- 
tals on a surface, a variety of patterns have been 
reported, ranging from ordered two-dimcnsional(2D), 
three-dimensional(3D) nanocrystal superlattices (NCSs) 
Gl (3 fractal- like aggregates 0, to more percolated 
networks 0, 0, ^| . A common notion is that the self- 
assembling process occurs at the liquid-substrate inter- 
face 01 1 where the competing effect from the diffusion 
of nanocrystals along the substrate and solvent dewetting 
could lead to pattern formation far away from equilib- 
rium |l2j |. Having such vast different macroscopic struc- 
tures of nanocrystal assembly hinders the exploration of 
their physical properties UjJ as well as creating 
obstacles for device applications Il5| . It is therefore 
essential to further understand and control the nanocrys- 
tal self-assembly mechanism. 

Recently, it has been shown that highly ordered 2D 
NCSs with domain size up to tens of microns, can be 
formed on a silicon nitride substrate by evaporating a 
nanocrystal colloidal droplet with slight excess amount of 
dodecanethiol ligand molecules in the solution [fj]. Fig- 
ure la shows a transmission electron microscopy (TEM) 
image of a region of such highly ordered NCSs. The high 
degree of ordering and large domain size indicates the 
formation of NCSs most likely occurs through nucleation 
and growth mechanism. It is difficult, however, to explain 
the high crystallinity of the superlattices through a mech- 
anism in which self-assembly occurs on the substrate at 
the last instance when the solvent dewets the surface. 



In this letter, we report an in situ, non-intrusive, small 
angle x-ray scattering (SAXS) measurement to elucidate 
the formation of NCSs as the colloidal droplet evapo- 
rates. We show that, contrary to the previous notion, 
these highly ordered 2D superlattices are formed at the 
liquid-air interface of the liquid droplet during the sol- 
vent evaporation. Changing the solvent evaporation rate 
can lead to either 2D or 3D NCS formation using the 
same colloidal nanocrystals. 

We used monodispersed gold nanocrystals with an av- 
erage diameter of 7.5 nm and 5.8 nm in two different 
experimental runs |l6l respectively, and with the 
particle number concentration adjusted to be approxi- 
mately 10 13 mL _1 , just enough to form a monolayer. The 
volume concentration of thiol was 0.63%. The SAXS ex- 
periments were conducted at the 1-BM beamline of the 
Advanced Photon Source (APS). 10 (A of colloidal so- 
lution was deposited onto a rectangular silicon nitride 
substrate (3 x 4mm 2 ). The incident synchrotron x-ray 
beam was monochromatized to 8.0 keV by a double- 
crystal monochromator while two sets of horizontal (H) 
and vertical (V) slits were used to define the beam size to 
0.2 (H) x0.2 (V)mm 2 . We chose a laboratory coordinate 
system as shown in Figure lb so that the incident x-ray 
beam is along the z-direction and the substrate surface 
is in the y-z plane. The samples were initially aligned 
using the substrate as a reference, thus the top surface 
of the substrate was centered in the x-ray beam along 
the x-direction. During the experiment, the substrate 
was slightly titled with an angle of 0.3° relative to the 
incident beam. The scattered x-rays passed through a 
helium flight path and were collected by either an image 
plate (IP) or a charge-coupled device (CCD) detector. 

The dynamical self-assembling process can be observed 
by positioning the incident beam at the bottom center of 
the droplet to accommodate the changing meniscus of the 
droplet as it evaporates. Figure 2 shows a typical time 
evolution of the SAXS patterns after the droplet is de- 
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FIG. 1: Characterization of 2D nanocrystal superlattices 
formed by colloidal droplet evaporation: (a) Transmission 
electron microscopy (TEM) image of a highly ordered gold 
nanocrystal monolayer. Inset shows Fourier transformation of 
the superlattice. (b) Schematic diagram of scattering geome- 
try used for the in situ SAXS measurement. The dimension 
of the substrate is exaggerated in order to clearly illustrate 
the experimental setup. 
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FIG. 2: (a-d)/n situ SAXS patterns of 2D NCSs formation 
during the droplet evaporation. Time in unit of minute is in 
reference to the deposition of colloid droplet. The average 
diameter of the nanocrystals is 7.5 nm determined by TEM. 
Cross bar feature in each frame is the shadow image of the 
beam stop, (e) Time-evolution of the normalized scattering 
intensity of the monolayer powder diffraction peak for (10) 
powder diffraction peak (solid line is a linear fit to the data) . 



posited on the substrate. The droplet is evaporated in air 
at room temperature (23°C) with the mass of the droplet 
decreasing approximately at 0.66mg per minute. The ini- 
tial thickness of the droplet is typically more than 2 mm. 
No visible scattering pattern is observed in the first two 
minutes indicating that, initially, there are no ordered 
superstructures in the droplet (Figure 2a). After that, 
an elliptically shaped diffraction ring is observed (Figure 
2b) . The section of the scattering ring near q y = grad- 
ually becomes more diffuse in the q^-direction and even- 
tually disappears, while the intensity of the scattering 
corresponding to the in-plane long-range order increases 
dramatically (Figures 2c-2d). The scattering pattern in 
Figure 2d remains unchanged for hours during which a 
thin liquid film, less than 100 /im in thickness, with a high 
concentration of dodecanethiol remains on the substrate 
surface. To illustrate the kinetics of the nanocrystal array 
growth, Figure 2e shows the experimental time evolution 
of the intensity of the (10) powder diffraction peak at 
q x = 0, plotted in a semi-natural log format. The scat- 
tering intensity increases exponentially with time before 
becoming saturated after 8 minutes (data not shown). 
When the liquid film is completely dried and the NCSs 
are deposited directly on the substrate, the downward 



scattering is blocked b y th e substrate in the grazing in- 
cidence geometry [lil Il9j. In contrast, the scattering 
patterns in the liquid film (Figure 2d) is highly symmet- 
ric in the vertical direction (q x ), indicating the 2D NCSs 
are formed at the liquid-air interface which is elevated 
above the substrate. The self-assembly of nanocrystals at 
the liquid-air interface was further confirmed by using a 
wider x-ray beam measured to be 1.5 (H)x 0.2 (V)mm 2 , 
purposely positioned near the edge of the droplet [2fj| . 
Bragg diffraction patterns from the 2D NCSs were ob- 
served rotating along with the direct reflection of menis- 
cus surface as the droplet flattens due to evaporation. We 
should emphasize that, unlike experiments carried out on 
the Langmuir trough [2l| , no aqueous subphase was used 
in our experiments. Both the nanocrystal ligand and the 
toluene solvent are hydrophobic in nature. 

A comprehensive understanding of the formation of the 
2D NCS monolayer, including the early stage scattering 
ring, can be achieved by comparing the simulated SAXS 
patterns during the entire evaporation process with the 
experimental data. The initial elliptically-shaped x-ray 
scattering patterns as well as its evolution can be mod- 
elled by the "powder" diffraction of 2D superlattice do- 
mains with varying x-ray incident angles. As shown in 
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the schematic cross section diagram in Figure 3A, the 
x-ray incident angle a at the liquid-air interface varies 
during the droplet evaporation (while the edges of the 
droplet are pinned at the substrate). To model the x- 
ray scattering patterns from the 2D lattices with their 
lattice plane orientation changing with time, we define 
a sample coordinate system (x',y') that is attached to 
the lattice plane as shown in Fi gure 3a. In a kinemati- 
cal approach of x-ray diffraction [22| , the averaged pow- 
der diffraction intensity from many 2D lattice domains 
in the reciprocal space is distributed in a circular ring 

defined by J q 2 x , + q 2 , = |qh,k|, where (q x >,q y i) arc the 

momentum transfer (MT) components in the sample co- 
ordinates which satisfies the Bragg diffraction condition 
defined by the scattering vector qh,k- The scattering 
patterns recorded by the 2D detectors in the laboratory 
space are related to the circular powder ring in the sam- 
ple space by a linear transformation between the two co- 
ordinates. Therefore the Bragg scattering condition in 
the laboratory coordinate is described approximately by 

y/ (q x sin a) 2 + q% = |qh,k|, where (q x ,q y ) are the MT 
components in laboratory coordinate. The short axis of 
the ellipse is in the g y -direction and has a magnitude of 
|qh.k| , whereas the long axis is in the ^-direction and has 
a magnitude of |qh,k|/sina . For example, the elliptical 
scattering pattern in Figure 2b indicates, at 3.5 minutes 
into the evaporation process, the 2D NCSs formed at the 
interface which is orientated about 52.4° with respect to 
the incoming x-ray beam. The complete scattering in- 
tensity is the product of the lattice structure factor and 
the form factor of a spherical nanocrystal expressed as 

\F(q)\ = | [ 3( S in( g r)-^co g ( 9 r)) ]|2| ^ g = J q 2 + q 2 + q 2 

and r is the radius of the sphere. As the liquid droplet 
evaporates, the scattering geometry approaches a grazing 
incidence angle with respect to the sample (a = 0.3°). In 
this case, the scattering intensity distribution along the 
^-direction corresponds to the structure factor of the in- 
plane lattice and that along q x reveals the form factor of 
the individual nanocrystal. Figures 3b-e show the simu- 
lation results at different stages of evaporation, based on 
a 2D hexagonal close-packed monolayer of nanocrystals 
with an average diameter of 7.5 nm (polydispersity of 
7%) and a 10.3 nm center-to-center spacing. Because the 
incident angle can be deduced from the simulation, it is 
possible to obtain the NCS scattering intensity normal- 
ized by the x-ray beam footprint size. This allows a pre- 
cise determination of the self-assembly kinetics. The x- 
ray scattering intensity data shown in Figure 2e is found 
to be proportional to the integrated intensity after being 
normalized by the changing x-ray footprint size on the 
sample due to the variation of x-ray incident angle. 

The formation of the 2D NCSs can be understood using 
a kinetic " crushing" model proposed by Nguyen and Wit- 
ten [2^] . The diffusion constant {D) of the nanocrystals 
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FIG. 3: Simulation results of the SAXS experiments, (a) 
Schematic cross sectional diagram of the x-ray beam and the 
colloidal droplet, (b)-(e) Simulation results based on a mono- 
layer of 2D domains oriented at an angle a with respect to 
the incident x-ray beam at the liquid-air interface. The time 
in each frame corresponds to the real time used in our exper- 
iments. 

with a hydrodynamic radius R (5.1 nm) can be estimated 
using the Stokes-Einstein relation (D = UbT /&TrqR) and 
the viscosity of dodecanethiol (77 = 2.98cp), assuming 
that toluene completely evaporates near the liquid-air in- 
terface. Since the major change in scattering patterns in 
Figure 2 occurs at a time scale of one minute, we can com- 
pare the diffusion distance of nanocrystals in one minute 
with the change of liquid-air interface during the same 
time interval. Performing a random walk, the average 
distance a nanocrystal can diffuse within one minute is 
V < r 2 > = y6Dt « 70/im . Since the air-liquid interface 
decreases faster than 70 /im/min during the initial evap- 
oration, as measured by monitoring the thickness change 
of the droplet in real-time, nanocrystals will accumulate 
at the 2D liquid-air interface and eventually induce 2D 
crystallization. 

The exponential increase of the scattering intensity 
(Figure 2e) also supports the interface " crushing" model, 
in contrast to the power law behavior in the conventional 
2D diffusion coarsening mechanism [24j]. In the "crush- 
ing" case, the domain growth mainly occur through in- 
corporating nanocrystals impinging the domain surface 
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FIG. 4: In situ scattering patterns obtained during a slow 
evaporation of the colloid droplet. In this experiment, the 
average diameter of nanocrystals is 5.8 nm determined by 
TEM. These patterns were recorded by the CCD detector. 



tal accumulation at the liquid-air interface, whereas un- 
der a much slower evaporation condition, nanocrystals 
can diffuse away from the interface and 3D NCSs can 
eventually form inside the droplet. The kinetics driven 
self-assembly process observed in our nanocrystal col- 
loidal system, might also play an important role in other 
complex systems. Exploring these systems using real- 
time x-ray scattering techniques should enable us to bet- 
ter control the organization of nanoscale building blocks, 
and also study many interesting phenomena, such as par- 
ticle and domain dynamics and order-to-disorder phase 
transition. 

We thank T. T. Nguyen, T. A. Witten, T. P. Bigioni, 
H. M. Jaeger and C. M. Sorensen for extensive discus- 
sions. This work was supported by the U.S. Department 
of Energy (DOE), BES-Materials Sciences, under Con- 
tract #W-31-109-ENG-38, by DOE Center for Nanoscale 
Materials, and by the University of Chicago - Argonnc 
National Laboratory Consortium for Nanoscience Re- 
search (CNR). The use of APS is supported by Office 
of Science of DOE. 



from below when it collapses with the evaporating liquid- 
air interface. The collected nanocrystals either fill the va- 
cant sites or move to the edge of the domain. The growth 
rate of a domain with size S is dS/dt^nSv , where n is 
density of nanocrystal and v is the collapsing rate of in- 
terface. We note approximately n and v remains constant 
during most of the evaporation process. The " crushing" 
effect thus leads to an exponential growth of domain size, 
which is a much more dominant process than the power- 
law growth by 2D diffusion at the interface. 

One prediction of this kinetic " crushing" model is that 
the formation of the 2D NCSs should be affected by the 
initial evaporation rate. A slower initial evaporation rate 
could allow the nanocrystals to diffuse away from the in- 
terface before the nanocrystal concentration reaches its 
critical value for 2D crystallization. The formation of su- 
perlattices will then occur in the bulk of the droplet once 
the overall concentration of nanocrystals in the droplet 
exceeds a critical concentration for 3D nucleation. Figure 
4 shows the time evolution of scattering patterns when 
the droplet is evaporated at approximately 0.22mg per 
minute, much slower than the ambient condition. In- 
deed, 3D supcrlattices form and persist during the entire 
evaporation process. Using the same colloidal solution, 
but through varying the evaporation rate, both the 2D 
NCSs and 3D NCSs formation can be reproduced consis- 
tently. 

Our measurements suggest the kinetics of evaporation 
can strongly affect the structure of NCS formation. 2D 
NCSs can form exponentially at the liquid-air interface 
when the evaporation is fast enough to induce nanocrys- 
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